Coulomb interacting Dirac fermions in disordered graphene 
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We study such experimentally relevant characteristics of the Coulomb interacting Dirac quasipar- 
ticles in disordered graphene as the quasiparticle width and density of states that can be probed 
by photoemission, magnetization and tunneling measurements. We find that an interplay between 
the unscreened Coulomb interactions and pseudo-relativistic quasiparticle kinematics can be best 
revealed in the ballistic regime, whereas in the diffusive limit the behavior is qualitatively similar to 
that of the ordinary 2DEG with parabolic dispersion. 



The recent advances in microfabrication of graphitic 
monolayers^ have made it possible to test and confirm the 
earlier theoretical predictions of anomalous, relativistic- 
like, kinematic properties of the electronic states in 
graphene^. Thus far, both experimental and theoreti- 
cal studies have been primarily focusing on the unusual 
(magneto)transport phenomena which, for their most 
part, can be described in terms of non-interacting Dirac 
quasiparticles propagating ballistically in the bulk or in 
finite geometries. 

However, it has long been recognized that a (nearly) 
degenerate two-dimensional semimetal such as graphene 
might provide a unique playground for studying the ef- 
fects of the Coulomb interactions, because the latter are 
expected to remain essentially unscreened^. 

Among the previously discussed manifestations of 
the Coulomb correlations is a possible opening of the 
interaction-induced excitonic gap at sufficiently strong 
Coulomb couplings^. Alternatively, such a gap can also 
be generated by a magnetic field, which phenomenon rep- 
resents the Dirac counterpart of FQHE in the conven- 
tional 2DEG with parabolic electron dispersion, as was 
pointed out in Ref.- (see also Ref.^ for the references to 
the earlier studies of a related phenomenon of " magnetic 
catalysis" in abstract field-theoretical setting). 

Should a spectral gap develop, it would be exhibited 
by the conductivity and other transport characteristics. 
However, apart from the recent observation-^ of a com- 
plete lifting of the four-fold degeneracy of the n — 
Landau level (which phenomenon would indeed be con- 
sistent with the scenario of a field-induced gap opening, 
resulting in a spontaneous breakdown of the sublattice 
symmetr y^i^) , no conclusive evidence of such a behavior 
has yet been found. 

Also, if proven to be of a genuine bulk nature (as op- 
posed to being due to magnetic impurities, edges and/or 
structural defects), the previously reported weak, albeit 
robust, ferromagnetism in pyrolytic graphite^ could be 
indicative of a possible instibility towards a (weakly) 
ferromagnetic excitonic state with unequal gaps for the 
spin-up and spin-down electrons'*. 

Obviously, a further experimental work is needed in 
order to accertain a real status of the scenario of a latent 
excitonic insulator proposed in Rcf.*, as well as contrast- 
ing it with such alternative predictions as that of the 
Stoner instability resulting in a fully polarized ferromag- 



netic state^.. 

In light of this uncertainty, in the present Com- 
munication we focus on the effects of the moderately 
strong Coulomb correlations which might not be powerful 
enough to generate a finite gap in the Dirac spectrum. As 
we demonstrate below, even in this case the quasiparticle 
properties can be affected in a number of experimentally 
relevant ways. To that end, we study the quasiparti- 
cle width and density of states in both, the ballistic and 
diffusive regimes, and contrast the results against those 
pertaining to the ordinary 2DEG with parabolic disper- 
sion. 

An extensive experience gained in the course of the 
previous studies of the conventional 2DEG suggests that 
such transport characteristics as the longitudinal and 
Hall DC conductivities may not provide the best means 
of revealing the Coulomb correlations. A general rea- 
son is that the two-particle response functions probed 
by transport measurements appear to be only weakly af- 
fected by such correlations due to a routine cancelation 
between the (potentially, large) fermion self-energy and 
vertex corrections. In that regard, a greater insight into 
the physics of interacting Dirac fermions can be provided 
by various single-particle probes, including photoemis- 
sion, tunneling, and magnetization measurements. 

The low-energy properties of graphene are governed 
by the electronic states in the vicinity of one of the two 
inequivalent nodal points (a = 1,2). Such states can be 
described by the Dirac Hamiltoniar^i^ 

F = iT;^ ^ /"*t^[<T,V, + (-l)"a^V,]*„ (1) 
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where vp is the Fermi velocity, go = 27re^/eo''^F ^ 3 is the 
bare value of the dimensionless Coulomb coupling, and 
(Ji is a triplet of the Pauli matrices acting in the space 
of (pseudo)spinors *„ = [tpa{A),'4'a{B)) composed of 
the values of the electron wave function on the A and B 
sublattices of the bipartite hexagonal lattice of graphene. 
The effects of the Coulomb interactions on the fermion 
propagator and interaction function are encoded in the 
fermion self-energy S and polarization operator H 

G'^(W, P)~^ = (W + ll)i - VFi&xPx + {-l)°'ayPy) 



+E«(l^,p), y«(c^,q) = [^+n«(^,q)]-i (2) 
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In the one-loop approximation, the former is given by the 
expression 

t%,p)= (3) 

— ^[/my4(c., q)G«(e + CO, p + q) coth — 

e + UJ. 



V^iuj, ci}ImG"'{e + t^, p + q) tanh ■ 
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The hnear-in-momentum term in ReY.^ {e, p) gives rise to 
a renormahzation of the running Coulomb coupling g{uj) 
described by the RG equation derived in Ref J. 
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d\n{VL/uj) 
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(4) 



The solution g{uj) « 5o/[l + (so/Stt) ln(r2/ijj)], where fi is 
an upper cutoff of order the electronic bandwidth, shows 
that the effective Coulomb coupling slowly decreases with 
decreasing energy. 

At a finite temperature T, chemical potential /x or elas- 
tic quasiparticle width 7 the RG flow described by Eq. (4) 
terminates below the energy scale ~ max\T, /i, 7]. 

At T > and/or in the presence of disorder, a func- 
tional form of the fermion polarization operator becomes 
quite prohibitive. However, in the ballistic limit and near 
half- filling (/i « 0), it can still be approximated as follows 
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), Q+«r (5) 
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where Q^ = c 

In Eq.(5), the first expression is the temporal compo- 
nent of the Lorentz-invariant free fermion polarization 
bubble computed at T = 0, whereas the second one 
exhibits the (well known in high energy physics) phe- 
nomenon of " thermal Debye screening" and concomitant 
Landau damping due to thermally excited quasiparticles. 

A straightforward analysis of Eq.(3) shows that 
the inelastic quasiparticle width defined as r(e,p) — 
ImTrY,f^{e, p) exhibits a strong "light-cone" singularity, 
akin to that previously encountered in the studies of the 
normal quasiparticles in d-wave superconductors where 
the commonly quoted T'^-behavior of the inverse quasi- 
particle lifetime represents a rough estimate that is only 
applicable to the thermal quasiparticles with energies and 
momenta e ^ vp ^ T (see RefJ^). 

When evaluated to the lowest order in the Coulomb 
coupling, r(e,p) appears to be discontinuous ai e — vfP 
and singular at e, vfP — >■ 0: 



STie,p) ^ g^eiPl)P+, P+>T 
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where P| = e^ ± VpP^, and d{x) is the Heaviside step- 
function. 

In both regimes, the dominant contribution comes 
from the transferred momenta of order ~ P-(_ due to 
the long-ranged nature of the Coulomb coupling (for a 
screened interaction, the last line in Eq.(6) would be re- 
placed with ~ (rVP+)i/2). 

At p = T = the linear energy dependence of the 
quasiparticle width was predicted in Refs.'^. It is worth 
noting, however, that in the case of bulk graphite dis- 
cussed in Refi^ the linear dependence would only hold 
at the momenta higher than the inverse inter-layer sep- 
aration 1/d, while for q '^ 1/d the screened interaction 
potental becomes less singular, ^(q) « goid/qY^'^. 

By contrast, in the case of graphene Eq.(6) would 
indeed hold all the way down to the energies ~ 
max[T, fi,"f], if it were not for the higher order correc- 
tions. In order to estimate their effect, we recalculate 
Eq.(3) with the RPA-dressed interaction function ac- 
counting for the fermion polarization (5), which proce- 
dure yields the result 

r(e,p)^0(p2)^lng, P+>gT 

^9iPl)ig^y/^ln^, T<P+<gT (7) 

while for P+ < T the logarithmic factor (which would 
have been absent altogether for any interaction less sin- 
gular than Coulomb) disappears from Eq. (7) . In contrast 
to Eq.(6), the result (7) remains continuous at the thresh- 
old e = VfP- 

In the presence of potential disorder, the low-energy 
quasiparticle width is dominated by the elastic part of 
the self-energy. A closed expression for the latter can 
be readily obtained in the case of short-range impurities 
with concentration n^ and scattering amplitude u 



S^(e,0) 






(8) 



Eq. (8) allows for a self-consistent calculation of the zero- 
energy quasiparticle width 7 = /rnTrE^ (0,0), ranging 
from the Born (u —f 0) to the unitarity (w ^ cxd) limit. 
It is worth mentioning, however, that Eq.(8) would need 
to be further modified in the potentially relevant case of 
the Coulomb impuritiesii. 

Having obtained 7, one can compute the non- 
interacting DOS at the Fermi energy and the correspond- 
ing Drude conductivity 



vci = ImTr 



P 



<(0,P) 



^27rw|, 



In 



o-Q 



' ~2 i 

7 vfn 
r4 Ml 



-—maxl , J 
h TT 7 



(9) 



where (5^(6, p) accounts for the impurity- induced broad- 
ening, but does not include any inelastic scattering. 



In order to study a crossover between the ballistic and 
diffusive regimes, we use the formula 



n^(^, q) 



1 



4wf ^/vj^^~~iuj'+lrf)^ - 7 



(10) 



interpolating between Eq.(5) for Q^ ^ 7 and the 
standard diffusive expressions, such as 11^ (a;, q) = 
CToq^/(-Dq^ — iuj), for Q+ < 7 (here D = CTo/t'o is the 
diffusion coefficient). 

In the ballistic limit (e, T ^ 7) the total self-energy 
is approximately given by the sum of Eqs.(7) and (8), 
whereas in the opposite, diffusive, regime (e, T ^ 7) the 
inelastic width can be found from a self-consistent equa- 
tion 



r(e,T)= / 
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whose solution behaves as 



r(e,r)^"^°"[^'^]in. ^"-"'^ 
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(11) 



(12) 



where g^= g{u} ^^). 

The total quasiparticle width can be deduced from the 
ARPES data^^. Alternatively, it can be inferred from the 
dHvA experiments and (to the extend that the oscillating 
part of the resistivity is indicative of the behavior of the 
single-particle Green function) the SdH ones. Namely, 
fitting the magnetization data to the formula 



AM{B) 



nsmjimn'^/eB) 2nn^,^i^,.T)/eB 



eBvj, ^_^ smh{2Tr'^nT fi/eBY 



(13) 

can provide, apart from such a spectacular hallmark of 
the free Dirac kinematics as the geometric (Berry) phase 
ttM,^ a valuable information on the energy/temperature 
dependence of the quasiparticle width. 

In the same spirit, one can estimate the Coulomb- 
controled phase breaking time T^{T) ~ {T/ao)lnaQ 
whose temperature dependence (familiar from the theory 
of the conventional 2DEG) can be manifested by magne- 
toresistance associated with the localization corrections 
to the Drude conductivityi^. 

Yet another viable experimental probe is provided by 
tunneling measurements. Previous studies have been pri- 
marily concerned with the behavior of the electronic DOS 
in the vicinity of strong potential impurities^^. However, 
despite ofTering a greater experimental observability of 
such prominent features as a resonant peak at (or close 
to, if the particle-hole symmetry is broken by subdomi- 
nant terms in Eq.(l)) zero energy, the near-impurity DOS 
appears to be highly non-universal and, therefore, reveals 
more information about the impurity potential itself than 
about the Coulomb correlations in the host electronic sys- 
tem. Notably, a typical plot of the near-impurity D0S1& 



appears to be very similar to that obtained in the case 
of a d-wave superconductor (see, e.g.r^ and references 
therein) where the Coulomb interactions would be com- 
pletely screened out by the condensate. 

In view of the above, in what follows we concentrate on 
the bulk DOS, the first interaction correction to which is 
given by the expression 

d,,{e) = -5]/mTr[Go^(e,p)]2E«(e,p) ^ 

77- f ^ 
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(14) 



where 7 = ^a'^g^. In the ballistic regime, the correction 
to the bare (linear) DOS features an additional (as com- 
pared to the case of the conventional disordered 2DEC ^^) 
logarithmic factor due to the aforementioned kinematic 
"light-cone" singularity, while in the diffusive limit one 
obtains the same diffusion-related (double-log) enhance- 
ment, as in the standard case. 

Associated with the DOS correction (14), there are the 
Altshuler-Aronov-type contributions to such observables 
as specific heat and quasiparticle conductivity which, un- 
like their weak-localization counterparts^, can not be 
readily suppressed by external in-plane magnetic field. 

Beyond the leading approximation, one finds an in- 
terference between the Coulomb interactions and disor- 
der, which further modifies the behavior of the idealized 
(clean and non-interacting) Dirac fermion system. Given 
the large bare strength of the Coulomb interaction, the 
higher order terms might contribute significantly, thus 
prompting one to employ an adequate non-perturbative 
technique. 

To that end, we make use of the tunneling action 
method of Refsi^. Adapting this approach to the case 
of graphene, we cast the tunneling DOS in the form 



Ke) 



-ImTr / G^(0,t)e-^W+"*di 
I" J-00 



(15) 



The disorder-averaged real-space/time Green function 



G^(0,t) 



oc e 
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jt is computed in the absence of the 



Coulomb interactions, while the latter are incorporated 
through the (imaginary part of) the action 

5W=/Jcoth^^/my(c.,q) 



X / dii / di2e"'"(*i"*"' < e*i(""(*i)-''(*^» > (16) 
/o "'0 

which describes the spreading of the excess charge asso- 
ciated with an act of tunneling into the graphene sample 
from, e.g., the STM tip. 

In the path-integral language, the averaging in Eq.(16) 
is carried out over all the quasiparticle trajectories r(i) 



contributing to the tunneling amplitude^^. In the balhs- 
tic regime, one obtains < e*''^''!^'"^) >« 1, whereas in the 
diffusive hmit < e'^^''^-''^) >« e-^'i\ 

To facihtatc a direct contact with experiment, we eval- 
uate the tunneling conductance 



G{V, T)-^J [niV + e) - n(e)]z.f ^(y + eHe)de 



de- 



v{e) 



'Tcos\i\V + e/2T) 



(17) 



where vpL{e) ~ const is the electron DOS of the normal 
(Fermi-liquid-like) STM tip biased at a voltage V . 
In the baUistic regime {V,T '^ 7) Eq.(16) yields 



S{t) 



_9L 



h\{nt), .go < 1 



^■ln(m)ln[— ln(m)], 50 > 1 



(18) 



TT^ e 

thereby resulting in the approximate power-law behavior 

G{V,T)cxmax[V,T]^+'^ (19) 

At weak coupling, the "zero-bias anomaly" (19) features 
a purely algebraic behavior with the anomalous exponent 
■q = gl/{4:n)'^. In contrast, at strong bare coupling the 
energy dependence of g{u!) gives rise to an approximate 
power-law decay where the effective exponent t](V, T) 
deviates slowly (only as ~ Inln^ /max\V,T\) from the 
universal value ry = (I/tt^) ln(87re) k, 0.43 attained at 
max[V,T] = n. 

In the diffusive regime {T,V < 7), the running cou- 
pling g{u}) levels off at the value g^, and the counterpart 
of Eq.(18) reads 



S{t) 



1 



(47r)2cro 



In(t7)[ln(t7) + 0(1)], I/7 < t < l/T, 
774^1n^[ln^ + 2rt], t>l/T (20) 



As a result, for /x ^ 7 (or ctq ^ 1) there exists an interval 
v^o ^ ln(7/maa;[y, T]) < ctq where one obtains the de- 
pendence similar to that of the conventional disordered 
2DEGi^ 



G(V^,T)cxi/oexp(- 
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' max\V, T] 

(21) 

At still lower biases and/or temperatures {m.ax\V,T] < 

^g-47r cro'j ^j^g conductance resumes a linear dependence 



)) 



G{V, T) ex 



g47r^CTo 



max\V, T] 



(22) 



reminiscent of the non-interacting DOS, but with a com- 
pletely different prefactor. In contrast to the interme- 
diate asymptotic regime (21) that can only occur in 
the strongly metallic case {jjl ^ 7), the linear depen- 
dence (22) might be expected to set in at the lowest bi- 
ases/temperatures for an arbitrary electron density. It 
is worth emphasizing, however, that, unlike in the case 
of the conventional 2DEG, the bare DOS of graphene is 
entirely due to disorder at low electron densities (/x < 7). 

In summary, we analyzed the effects of moderately 
strong Coulomb interactions on the Dirac quasiparticle 
excitations in graphene. Taken at their face values, the 
above results for the quasiparticle width and DOS sug- 
gest that the Dirac physics can be best revealed in the 
ballistic regime (see Eqs.(7) and (19)), while the diffusive 
dynamics of this system (see Eqs.(12) and (21,22)) ap- 
pears to be deceptively similar to that of the conventional 
2DEG. These predictions of both, novel and mundane, 
features can be tested in future experiments on photoe- 
mission, tunneling, and magnetization measurements. 
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